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I. MOTIVATION 



The possibility to create strangelets or droplets of metastable cold strange quark mat- 
ter in ultra-relativistic collisions has been proposed and studied by several authors 
The existence of such exotic states, as well as metastable exotic multi-strange baryonic 
objects (memo's), has fundamental importance in cosmological models and in the un- 
derlying description of strong interactions. Several experiments are being carried out at 
the Brookhaven AGS (E864, E878) and the CERN SPS (NA52), which look for strangelet 
production (searching for small Zj A ratios) [P|-p!U|. 

Among the different theoretical approaches a mechanism of separation of 

strangeness from anti-strangeness (distillation process) has been proposed @] during 
hadronization of a system at finite baryon densities. This scenario, which assumes a first 
order phase transition, predicts a relative time delay between the production of strange and 
anti-strange particles. In the present paper, we propose to search for this delay with the 



help of the novel correlation method proposed in [0,|T^]. It exploits the sensitivity of the 
directional dependence of the correlation function of two non-identical particles to the time 
order of their emission. This new technique can allow for the study of the transient strange 
quark matter state even if it decays on strong interaction time scales |p!3|-p!5|]. 

A different method, exploiting K^K^ correlations, also allowing for a study of the relative 
delays (but not their sign) in K and K emission was proposed in JlBl. It is based on the 
fact that the K^K^ system, due to the positive CP parity of K^, triggers out the symmetric 
combination of the K^K^ states, thus leading to the familiar Bose-Einstein correlation 
pattern [0. This pattern, contrary to the case of non interacting identical bosons, is 
however substantially modified by the effect of the strong final state interactions [|18|. 



Here we will use the general method of [pJ] , p!2| to study the impact of the time delays on 
the correlations in K^K~ pairs. The choice of this system, similar to the case of neutral 
kaons, apart from the possibility to determine the signs of the time delays, has the advantage 
of weaker distortions caused by the resonance production (more than 50% of kaons are pre- 
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dieted to originate from direct emission ^^), as compared to pions, while the experimental 
feasibility is better than in the case of neutral kaons. 

II. THE MIXED PHASE THERMODYNAMICAL APPROACH 

The dynamical evolution of the mixed phase consisting of a quark gluon plasma and 
hadronic gas will be described in a two-phase model which takes into account equilibrium 
as well as non-equilibrium features 

Within this model, two main assumptions are made. Firstly, the QGP is surrounded by 
a layer of hadron gas and equilibrium ( described by Gibbs conditions) is maintained during 
the evolution. Secondly, non-equilibrium evaporation is incorporated by a time dependent 
emission of hadrons from the surface of the hadronic fireball. Within this model, it is possible 
to follow the evolution of mass, entropy and strangeness fraction fg of the system and extract 
relative yields of particles which compare reasonably well with experimental data [|T^. 

One specific feature of these calculations is the prediction that the system quickly enters 
the strangeness sector leaving the usual fiq — T plane of the phase diagram. While the 
net-baryon number Ab decreases during the hadronization process, the s quark chemical 
potential /ig increases from /ig = MeV (for = if we start with a pure u,d quark 
phase) to several tens of MeV. Consequently, the strangeness fraction = {Ng — Ns)/Atot 
increases. This is the so-called "distillation process" which can result in the formation 
of stable or metastable blobs of strange matter in the case of low bag model constants 
5^/^ < 180 MeV. (In this cooling of the system is predicted.) Strange and anti- 

strange quarks, produced in equal amount in the hot plasma, do not hadronize at the same 
time. Because of the excess of massless (u and d) quarks (as compared to their anti-quarks) 
in the case of a baryon rich plasma, the hadronization of and (containing s quarks) 
will be favored preferentially to K~ and Kq ( containing s quarks). The evaporation from 
the surface of the hadron gas, which is rich of anti-strangeness, carries away K^, K^,... 
thus charging up the remaining mixed system with net positive strangeness. Fig.l shows, 
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for two rather different bag constants, the time evolution [|T9| of the baryon number A, the 
strangeness fraction fg and the temperature T for initial entropy S/A^^i = 10, 45 and fs = 0. 
Another prediction of interest is the hadronic freeze-out time. Strange and antistrange 
hadrons are not emitted at the same time since the hadron densities in the outer layer 
are dictated by the strongly time dependent chemical potentials and the time dependent 
temperature ||20|JT9|1 . One should also notice that the difference between N and K~ N 
cross sections in a baryon rich gas will slow down the diffusion of negative strangeness (s). 



thus leading to an additional creation of time differences P0| , p!9 . 

In the present work, this time separation characterizes the transient existence of the 
plasma. In the following, we will show that it can be used to search for the distillation 



process by using correlation techniques of unlike particle pairs 1 22, 23, 18 



III. UNEQUAL PARTICLE CORRELATIONS: TIME-ORDERING SENSITIVITY 

Particle correlations at small relative momentum are mainly driven by strong and 
Coulomb final state interactions (in the case of unlike particle pairs). 

At high bombarding energies, this effect has been often treated with a well known size- 
independent Coulomb correction factor. Contrary, at low energies, the Coulomb effect be- 
comes the main tool for the study of this evolution due to larger time sequence of emission. 

Unlike particle correlations have been used experimentally and theoretically for more 
than 15 years [p^p2| , particularly for the study of heavy ion collision mechanisms around 
the Fermi energy p2|-|25[|. Particularly, evidence for three-body Coulomb effects on two- 
particle correlations was observed in |2^ and taken into account to describe a variety of 



two-particle correlation patterns using classical trajectory calculations |2^ or a complete 
three-body quantum approach ||Tl] , |26| in the adiabatic limit. The sensitivity of unlike particle 
correlations to the order of particle emission was pointed out and studied using both quantum 



TT|,[T2| and classical trajectory [^] approaches. In the framework of the former one and under 



the assumption of sufficiently small density in phase space, the two-particle correlation 
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function at a given relative c.m.s. momentum q = 2k is determined by the modulus squared 
of the two-particle amplitude averaged over the relative c.m.s. coordinates f* of the emission 
points. The sensitivity of the correlation function to the time delays, or generally to the 
space-time asymmetries in particle production, appears due to the dependence of the two- 
particle amplitude on the scalar product kf*. In particular, in the limit of large relative 
emission times t = ti — t2, v\t\ ^ r, the Lorentz transformation from the source rest frame 
to the two-particle c.m.s {v and 7 are the pair velocity and Lorentz factor): 

r*L = 7{rL-vt), r*j, = rT, (1) 

indicates that the vector r* —'yvt is only slightly modified by averaging over the spatial 
location of the emission points in the rest frame of the source. So the vector f* is nearly 
parallel or antiparallel to the pair velocity v, depending on the sign of the time difference t. 
The dependence of the amplitude on k ■ f* is thus transformed into the dependence of the 
correlation function on — A; ■ v'~ft. Therefore, the sensitivity to the sign of the time difference 
t is due to the odd part, in k ■ v, of the correlation function. The mean relative emission 
time, including its sign, can be determined by comparing the correlation functions R'^ and 
R~ corresponding to A; ■ -u > and < |12|. Noting that for particles of equal masses 



the sign of k ■ v coincides with the sign of the velocity difference Vi — V2, we can see the 
simple classical meaning of the above selection. It corresponds to the intuitive expectation 
of different particle interaction in the case when the faster particle is emitted earlier as 
compared to the case of its later emission. 

For charged particles characterized by a large Bohr radius, \a\ 3> (r*), |Re/| (/ is the 
amplitude due to the strong interaction), the ratio (1 -|- R~^)/{1 + R~) at very small q takes 
on a simple analytical form p8|: 



(1 + i?+)/(l + R-)^1 + 2{rl)/a 1 - 2{-fvt)/a, (2) 

where the arrow indicates the limit v\t\ ^ r. Thus, for K'^K^-sjstem (a = —111 fm), each 
fm in the asymmetry (rj) or (7ft) transforms in a 2% change of the (1 -|- R'^) /(I + R~) ratio 



at g ^ 0. The correlation functions R'^ and R~ as predicted in |[TT|JT^ well agree with the 
proton-deuteron correlation functions measured in low energy nuclear collisions studied at 
GANIL p9| , p^ , pO| , |3T| . As proposed in |]12| , this technique can be extended for a study of the 
emission time differences between any kinds of interacting unlike particles. 

In preliminary studies using simple event generators, it has been demonstrated [|I^-[T^ 
that K^K^ pairs with mean emission time differences as small as ±5fm/c can give rise to 
observable differences between R'^ and R~ correlation functions. It has been also found that 
most of the effect at g < lOMeV/c originates from the Coulomb interaction between the 
two kaons and that, in full correspondence with Eq. (2), the later (on the average) emission 
of K~'s is associated with the (1 + R^)/{1 + R~) ratio less than unity. 

In this paper, we will use this method in conjunction with the dynamical two phase 
description (see previous section) in order to quantitatively demonstrate its sensitivity to the 
predicted delays between the and K~ emission related to the production of a transient 
strange quark matter state. 



Three sets of parameters [T^JH] have been chosen: 
set 1 ^ 2 (representing Au+Au at AGS energies [|I^): 

initial mass A = 394, entropy per baryon S/A = 10, initial net strangeness fs = 0; 

set 1: B^/^ = 160 MeV, strangelets (albeit metastable) are formed [|19|] . 

set 2: B^l^ = 235 MeV, strangelets are not distilled |jl9|- They are not (meta-)stable. 



set 3 ( representing S+Au at SPS energies |T9 



bag model constant B^^^ = 235 MeV, initial mass A = 100, entropy per baryon S/A = 45, 
initial net strangeness /s = 0. 



In the case of the low bag constant (B^^^ = 160 MeV) - set 1 - a cooling of the system 
is predicted leading to rather long kaon emission times. In the early stage mainly K~^^s are 
emitted (see Fig. 1, l.h.s. column) so that a cold strangelet emerges in a few tens fm/c. 
According to Eqs. (1), (2) and the corresponding discussion, the later emission of K~^s 
should lead to the correlation function ratio (1 + -R"'")/(l + R~) less than unity provided 
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that the asymmetry (r^) in K^K~ c.m.s. is dominated by the time asymmetry term (vt). 
As seen from table 1 and Fig. 2 this is indeed the case: the spatial asymmetry term {ri) 
contributes by less than 10% (in the same direction); the predicted value of the deviation 
of the correlation function ratio from unity of —26% at g ^ is in good agreement with 
the intercept value of the calculated ratio curve with the ordinate axis. Also, the deviation 
from unity remains important and statistically meaningful up to the region q 30MeV/c. 
Thus, in spite of experimental difficulties to measure correlations down to a few MeV/c, the 
signal proposed for consideration is accessible to usual correlation measurements. Since for 
various two-particle systems produced in the mid-rapidity region, the ordinary dynamical 
mechanisms lead to the intercept values deviating from unity by less than 10% [Q, we can 
consider the eventual observation of a large negative asymmetry in the K^K^ correlation 
function ratio as a signal of the strangelet formation. 

In the case of the high bag constant (B^/^ = 235 MeV) - set 2 and 3 - the system heats 
up slightly. This results in a fast hadronization and, subsequently, in a small difference 
between the and K~ emission times ||T^ (see the middle and r.h.s. columns in Fig. 1). 
Similar to the previous case the i^+'s are emitted earlier than i^~'s. Also the asymmetry 
(r2) is dominated by the time asymmetry term {vt)^ the spatial one, (r/,), contributing by 
about 25% in the same direction (see table 1, rows 2 and 3). However, the asymmetry effect 
is now much weaker (Fig.2): the predicted intercepts with the ordinate axis of the correlation 
function ratio deviate from unity by —4.7% and —2.3% for sets 2 and 3, respectively. 

Clearly, in this case, such a relatively weak asymmetry effect cannot be considered as a 
signal of the strangeness distillation without detailed studies of the impact of the ordinary 
mechanisms (rescattering effects, for example). 

Regarding the width of the effect seen in the correlation function ratios (Fig. 2), it 
decreases with the increasing (r*) (see table 1) in correspondence with the narrowing of 
the correlation effect itself. Thus, for set 1 ((r*) 18 fm) the effect rapidly vanishes with 
increasing q and extends up to about 30MeV/c only. For sets 2 and 3 {{r*) 6fm in both 
cases) the deviation from unity is softer in the same range of q values and extends up to 
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about 60MeV/c. 



IV. CONCLUSIONS 

We have presented a novel method which can be apphed to characterize the possible ex- 
istence of a strange quark matter distillation process in ultra-relativistic heavy-ion collisions. 
The method is based on the predictions [ll8| , |rTHl^] for unequal particle correlations and ex- 



ploits the predicted properties of the transient, strange quark matter state 0,^,^|T9 



even 



if it decays on strong interaction time scales. Using the description of this strangeness distil- 
lation process by a dynamical evolution model for the mixed QGP-hadronic phase, we have 
quantitatively demonstrated the sensitivity to the bag constant correlated to the stability 
of the quark matter possibly encountered in available experimental situations. 

For the lower bag constant case a strong and sharp negative asymmetry effect is predicted 
in the correlation function ratio, thus offering a clear signal of the predicted strangeness dis- 
tillation process provided the experimental resolution and statistics are sufficient to measure 
the ratio down to about 10 MeV/c. Weaker though wider signals are predicted for the higher 
bag constant hypothesis. In this case the competition with other collision mechanisms gener- 
ating the asymmetry is possible. These mechanisms, including the expansion of the hadronic 
system, are currently under investigation [32,33]. 

Finally, we think that the unlike particle correlation technique may offer a very valuable 
tool to disentangle between different and presently debated scenarios for the phase transition 
in QCD matter at finite net baryon density. 
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set 


{vt) (fm) 


{tl) (fm) 


{rl) (fm) 


(r*) (fm) 


1 


-10.9 


1.1 


14.3 


17.5 


2 


-1.4 


0.5 


2.6 


5.5 


3 


-0.7 


0.2 


1.3 


6.3 



TABLE I. Mean values of the relative space-time coordinates for pairs with relative momenta 
q < 50MeV/c calculated within the mixed phase thermodynamical model (see text). 
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FIG. 1. Time evolution of baryon number A, strangeness fraction fg, temperature T and yields 
of if+ and K" for set 1 (left column), set 2 (central column) and set 3 (right column) (see text). 

FIG. 2. Ratios of correlation functions (1 + i?+)/(l + R") for three parameter sets in a 
mixed-phase thermodynamical model (see text). The strange quark matter distillation process 
predicted in the case of a low bag constant (upper curve) results in a measurable deviation from 
unity in the range q < 30MeV/c. 
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